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CN . Abstract 

>>: 

^ ' The method of QCD sum rules is used to calculate the mass of a hypothetical stable dihyperon H with 

the quantum numbers of two A hyperons, whose existence was predicted by R. L. Jaffe [Phys. Rev. 
if) [ Lett. 31, 195, 617(E) (1977)] in the framework of the MIT quark-bag model. Within the accuracy of 

the method of QCD sum rules, which in the present case is ~ 20%, the results obtained here agree 
with those of Jaffe. However, the method of sum rules does not make it possible to determine whether 
'■ the mass of the dihyperon H lies above or below the AA threshold (i.e., whether H is stable). 
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The theoretical investigation of multiquark states {Q^Q^-, n + m > 3) and the experimental search 
Q^' for them may provide important information about the properties of the interaction of quarks at large 

distances. In particular, major efforts, both theoretical 0121 and experimental ^ have been directed 
D \ at the study of the dibaryon spectrum (n = 6, m = 0). 

' In this connection, interest attaches to Jaffe's prediction ^ that there exists a stable six-quark 

s-wave state - a dihyperon which is a singlet with respect to both color and flavor (with strangeness 
^ \ —2) and has = O"*" and mass 2150 MeV. The quantum numbers of H are identical to the quantum 

$H ' numbers of AA [two A(1115) hyperons], and its mass is smaller than the sum of the masses of two 

A hyperons. Therefore H can decay only through the weak channel (consequently, it is stable with 
respect to the strong interactions). 

The prediction that a dihyperon exists was obtained in the framework of the MIT quark-bag 
model [1]. It is important to test this model prediction solely on the basis of the fundamental principles 
of QCD. Such a test can be made by means of the method of QCD sum rules, using either the technique 
of finite-energy sum rules [H] or that of Borel sum rules [HI. The method of QCD sum rules is based 
on the fundamental principles of QCD, and it has proved its effectiveness in calculations of the masses 
of mesons (Si tZil iSl iSj and baryons [HI El ^ • 

In the present paper, assuming the existence of the above-mentioned dihyperon, we calculate its 
mass (with the aim of testing the prediction of Ref. 0) in the framework of the Borel version of the 
method of QCD sum rules 1^. 

According to the method of QCD sum rules, to calculate the H mass we must consider the 
correlation function 

n(g) = i J d^xe'''''{0\T{h{x)h+{0)\0). (1) 
Here h{x) is a scalar local current with the quantum numbers of the dihyperon H, constructed from 
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six quark fields, such that its projection onto the real dihyperon state \H{p)) (we are working with 
the assumption that this state exists) is nonzero: 



{0\h{0)\H{p)) = Xh, p^ = ml 



(2) 



Since such a current h{x) is not unique, the question of its optimal choice immediately arises. We recall 
that the problem of choosing the current already arose jlUUllj in the case of baryons, when the currents 
are constructed from three quark fields. In our case, when the current is constructed instead from six 
quark fields, this problem is much more complicated, since the number of independent currents with 
the given quantum numbers is much greater. The treatment of the current h{x) in the most general 
form, i.e., in the form of a linear combination of all the independent local operators with the quantum 
numbers of i/, is a very cumbersome problem. Therefore we confine ourselves to the choice of a few of 
the simplest currents with the required quantum numbers and analyze the dependence of our results 
on the properties of the studied currents. As in the case of mesons and baryons, we shall construct 
the current h{x) from quark fields without derivatives. 

We shall consider only the following six simplest (from the point of view of the calculation of the 
correlation function 11) currents with the quantum numbers of the dihyperon H (we recall that H is 
a singlet with respect to both color and flavor): 



/ll(x) =e^i^2^^e^4^5^«ea,a,a3e„,„3„,(VlCV'2)(V'3C75V'4)(V'5CV'6) , 

h2{x) = e^^^2^^e^^^^^«eaia2a3£a4a5a6(^lC75V'2)(V'3C75V'4)(^5C75V'6) , 
h3{x) = e^^^^^^e^*^^^«eaia2a3ea4a5a6(V'lC7^75^2)(V'3C75V'4)(V'5C7M75V'6) , 
hi{x) = e^^^2^*e^^^^^«eaia2a3ea4a5a6(V'lC75V'2)(V'3C75V'4)(V'5C75V'6) , 
h{x) = e^l^^^*e^^^^^«eaia2a3ea4a5a6(V'lC75V'2)(V'3C75V'4)(V'5C75V'6) , 



where ipn = V'a" (^) ^ quark fleld, a„ is the color index. An is the flavor index, the spinor indices are 
omitted, and a summation over repeated indices is understood; g^i^z^a a,nd £010203 are antisymmetric 
tensors which are invariant with respect to the flavor group SU-^ (3) and the color group SUc{2), 
respectively; C are the matrices of charge conjugation. [The choice of the current hQ{x) will be 
discussed below.] 

In accordance with the method of QCD sum rules, we shall calculate the correlation functions 



by means of Wilson's operator expansion, assuming that the vacuum expectation values of the local 
operators (the so-called condensates) are nonzero. The calculations must be performed in the region 
-g^ > 1 GeV^. In this region, on the one hand, the effective strong interaction constant ag is already 
small enough to calculate the coefficient functions of the operator expansion by perturbation theory; 
on the other hand, the nonperturbative corrections due to the nonzero vacuum expectation values of 
the local operators are already quite important. 

Because of the great complexity of the calculations of Ilj , we shall confine ourselves to a calculation 
in the zeroth order in a^. The fact that this approximation may be sufficient for calculations of the 
physical quantities in the framework of the method used here is confirmed by the practice of using the 
method, in particular by the calculation of baryon masses. We calculate Ilj in the linear approximation 



/l6(x) = /l5(x)-(l/3)/l2(x) 



(3) 




(4) 
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in the strange quark mass and in the zeroth order in the masses of the u and d quarks. The results 
of our calculations are as follows: 



III = 112 = T7:^3 = const • 114 
lb 



(5) 
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Here a = \167r'^ {qq)\, isotopic invariance gives {uu) = {dd) = {qq), and we also assume {ss) = {qq) for 
the strange quark condensate, since our analysis has shown that allowance for the difference between 
(ss) and {qq) practically does not alter our final results. In the expressions (5) we have retained only 
the terms which survive after application of a Borel transformation (see below); in the calculation 
of the vacuum expectation values of multiquark operators, we have made use of the hypothesis of 
factorization jlj; the factor l/(167r2)^ corresponds to the fact that the perturbative contribution is 
represented by five- loop diagrams. 

For each correlation function Ilj we can write the dispersion relation (below, we shall drop the 
index j, since the arguments are valid for all j) 

n(Q2) = r _ subtractions , = -q^ , (6) 

Jo s ~\~ 

where the spectral density p{s) = (l/7r)Imn(— s— iO) behaves as for large s, so that eight subtractions 
must be made in the relation (6). Applying a Borel transformation (i.e., the differential operator 



B = lim 



(n - 1)! V dQ 



QVn=y 



which makes the subtractions vanish) to both sides of Eq. (6), we obtain a Borel sum rule of the well 
known form 

U^{y) = - / e-'/yp{s)ds, (7) 

y Jo 
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where Il^{y) can be obtained in our approximation from (5) by the action of the operator B: 

BQ-.„g)=(-l)-.!/, = (8) 

In order to calculate the mass of the dihyperon, it is necessary to choose some phenomenological 
ansatz for the spectral density p{s). The simplest procedure which has proved to be successful in 
the calculation of meson and baryon masses is to represent p as the sum of the pole contribution of 
the unknown lowest state of the given channel and an effective continuum beginning at some effective 
threshold W and approximating the contributions of all the higher states. We shall also take a spectral 
density in the "pole plus continuum" form 

p(s) = Xj^5is - m\) + e{s - W^)p^'"'\s) , p^'"'\s) = (l/7r)Imn(-s - iO) (9) 

where the spectral density of the continuum is determined as the imaginary part of the complete 
approximation calculated for 11. Substituting the spectral density Q into Eq. Q and transferring the 
continuum contribution to the left-hand side, we obtain a sum rule for the correlation function Hi, 

^7 2 2„,4 1 4„, I ^ f 4 5 4 3 2 ^5 ^ 



-y a y a y + m, ay a y H a 

35^ 15 ^ 18 ^ V 9 ^ 27 ^ 324 y 

- - r e-'/yp'"''\s)ds = -X^e-"'H/y (10) 

y Jw^ y 

and a sum rule for the correlation function Xlg, 

43 7 1081 2 4 1081 4 /3034 . 2162 . o 209 .1\ 

— y' H a y H a y + ay^ H a:^y H a^- 

63^ 135 ^ 324 ^ V 405 ^ 243 ^ 5832 yj 

1 1 

- - / e-^V°"*(s)ds = -}?e-'^H/y (11) 

y Jw^ y 



where 



64 ' " 64 

The sum rules for 112,3,4 are identical to (10) apart from an unessential overall factor. We have 
also not written down the sum rule for the correlation function lis, since the coefficients in are 
numerically close to the coefficients in the correlation function [see Eq. (5)], so that the results 
obtained from the sum rule (11) for Hg are practically identical to the results of 

The sum rules must be considered in an admissible interval of values of y, in which the contribu- 
tions of the higher power corrections and, at the same time, the contribution of the effective continuum 
must not be too large. We shall call this interval $7 and impose the following restrictions for its deter- 
mination: 1) the continuum contribution is less than 50% of the calculated approximation for Il^{y) 
which means that the continuum contribution is smaller than the contribution of the dihyperon H to 
the sum rule; 2) the sum of all the nonperturbative corrections, excluding the leading nonperturbative 
term ~ (gg)^y^, is less than 50% of the calculated approximation for Il^{y). (Although with respect 
to the dimensions of y the term ~ ms{qq)'^y^ is the leading nonperturbative term, it is natural to 
include it in the constrained sum, since it is ~ rrig and numerically the case = should not differ 
fundamentally from the case of broken SU^ (3) symmetry, 7^ 0; moreover, if the term ~ ms{qq)'^y^ 
is not included in the constrained sum, this has little influence on the results.) 
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The choice of some a priori Umit (in our case, 50%) is highly conditional. In the final anaysis, its 
correctness is substantiated a posteriori by the consistency of the sum rule in the interval 0, which is 
found. 

Before particularizing the procedure of fitting the sum rules, we shall discuss the choice of the 
numerical value of the vacuum expectation value {qq) , which depends on the normalization point /Li. 
We shall choose /i^ ~ y G 17. Then the renormalization-group factor [as{y) /asin'^)]^" in the coefficient 
function of a given operator 0„ (where 5„ is determined by the anomalous dimenions of 0„ and of the 
current h) becomes ~ 1 for y € i}. Therefore we shall neglect the effects of the anomalous dimensions. 
In our case, the analysis of the interval shows that it is necessary to take ~ 1 GeV. If we begin 
with the value from Ref. 6, 

K^g)U~200 Mev = (240 MeYf = 0.014 GeV^' 

then the application of the renormalization group gives the first estimate \ {qq)\fir^i GcV = 0.024 GeV^. 
On the other hand, if we recalculate the value of {qq) according to the PCAC formula {qq) = 
— (l/2)/^m^/(mu + md), then by using the value m„ + m^l^^i ceV = 16 MeV which is now adopted 
for the masses of the light quarks ^2] we obtain the second estimate \{qq)\iir^i GeV = 0.011 GeV^. It 
can be seen that those two values differ by a factor ~ 2. 

Our final results depend little on the variation of {qq) in the indicated interval. This is perfectly 
natural. The value of {qq) is the only number which is built into our sum rule (inclusion of the terms 
~ nis should not substantially alter the results). Therefore it follows from dimensional arguments 
that the mass of the dihyperon is mn ~ {\{qq)\)^^^- Consequently, mn should not change by more 
than ~ 30% if {qq) varies by a factor 2. 

We must determine three parameters (m//, W, A) by fitting the sum rule in expression (10) [or (11)], 
i.e. on the basis of agreement between the left and right hand sides [which will be denoted by L{W, y) 
and R{mH,y), respectively] of the sum rule (10) [or (11)] in the admissible interval 0. We rewrite the 
sum rule in the form 

L{W,y)ye'^H/y = X\ (12) 

Let F{mH,W,y) denote the left hand side of Eq. (12). Since must not depend on y, we shall 
determine the values of mn and W by requiring the minimal dependence of F on y in the interval fl. 
Thus, only two parameters mn and W are effectively fitted, and A^ is determined simply as the mean 
value F{mH, W) of the function F in 0, for the values of mn and W which are found. We introduce 
a criterion for the determination of mn and W. We consider the relative deviation of F{mH ,W,y) 
from its mean value: 

F{mH,W,y)-F{mH,W) 
F{mH,W) 

We shall fix the values of mn and W by requiring the largest possible interval (the upper limit of O 
depends on W) subject to the condition that the relative error A does not exceed a few percent (i.e., 
the function F is practically a constant in fi, since it is a smooth function of y). For definiteness, we 
require that A < 5%. 

We analyzed the accuracy of the determination of the parameters on the basis of the fitting 
procedure described above as follows. The criterion for A was varied by a factor 2 in each direction 
from the given value A = 5%, and the value of the quark condensate {qq) was varied between the 
limits described above. The corresponding changes in the values of mn and W determined from the fit 
were not more than 20%. Therefore we can say that our sum rule fixes mn and W with 20% accuracy. 
The parameter A^ is determined with much lower accuracy with an accuracy up to a factor 2. The 



FL{W,y)-R{mH,y) 
R(mH,y) 



(13) 
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greater accuracy in the determination of mn and W than in that of A is natural, since mn and W 
appear in the arguments of exponential functions and, consequently, the sum rule is more sensitive to 
their variation. 

Our analysis has shown that the sum rule (10) for the currents hi — is poorly fitted (in other 
words, it is not possible to satisfy the criterion introduced above). At the same time, the sum rule (11) 
for the current /ig is well fitted over a rather large interval 0. In particular, for [qq) = —(240 MeV)^ 
the results of fitting the sum rule (11) are as follows. For = we obtain 

ruH = 2.0 GeV, W = 2.8 GeV, = 70 GeV^ , (14) 

and for = 0.2 GeV (Ref. 12) 

mu = 2.4 GeV, W = 3.2 GeV, = 150 GeV^ . (15) 

The y dependences of the left and right hand sides of the sum rule (11) for the parameter values (14) 
and (15) which we have found are shown in Figs, la and lb, respectively. 




Figure 1: Dependence of the left hand side (continuous curves) and right hand side (broken curves) 
of the sum rule (11) on y for the parameter values in (14) and (15): a) = 0; b) nig = 0.2 GeV. 

We recall the values obtained by Jaffe fl| for the mass of the dihyperon H (the values of the 
dibaryon masses obtained in the quark-bag model correspond to poles of the P-matrix jl3j and the 
physical states have masses 100-200 MeV lower than those predicted by this model): 

mH\ms=o = 1-76 GeV, m^lm^^o = 2.150 GeV . (16) 

It can be seen that the values (14,15) obtained for the dihyperon mass by the method of QCD sum 
rules agree within the accuracy of this method (in our case, the mass is determined with an accuracy 
(~ 20%) with the results ^ of the MIT quark-bag model. 

However, there immediately arises the question of what makes the current /ig, for which the 
results (14,15) are obtained (the results for /15 are practically identical to them), preferred over the 
currents hi — /14 [see Eq. (3)] which lead to the "poor" sum rule (10). 

In our opinion, the answer is that the current h^ satisfies the following two "physical" requirements. 
First, there exists a nonzero nonrelativistic limit for it [i.e. if the quark field ^{x) is represented in 
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the standard manner in terms of the small and large components, the term containing only the 
large component will be nonzero]. On the other hand, the currents hi and /14 do not possess a 
nonrelativistic limit. We note that in Refs. ^01 and ^Ij it was already pointed out that the existence 
of a nonrelativistic limit is a desirable property for the construction of currents in working with the 
method of QCD sum rules. 

Second, we represent the current /ig in the form of a product "singletiS'singlet" with respect to 
color, i.e., in the form = ^'i(x)r\I'2(x), where each of the operators ^1, and ^2 is a color singlet 
and r is some combination of Dirac 7 matrices. Then we obtain the representation (the currents hj 
are as a whole singlets with respect to both color and flavor) 

heix) = {^il{x)riC^,t{^il{x))0 (17) 

where the color-singlet operator ^' is a flavor octet: 

i^it)- = ^a,a,aAriC^^ri){rB\).e^'^''''{^%^^^ - (l/3)5g<^) . (18) 

Here a and [3 are spinor indices, and the remaining notation is the same as in (3). 

We note that the flavor octet (|18|) has the quantum numbers of the baryon octet and has been 
used to calculate the characteristics of baryons by the method of QCD sum rules ^3^5- Thus, if 
the current is represented in the form "singlet^singlet" with respect to color, it has the structure 
"octet®octet" with respect to flavor (each color singlet is a flavor octet). Since the hyperon A is a 
member of the baryon flavor octet and the quantum numbers of the dihyperon H are identical to the 
quantum numbers of AA, it seems natural to construct a flavor singlet h{x) in the form "octetCgioctet" . 

Physically, we can imagine the following picture. If we split the colorless dihyperon H into two 
colorless clusters and separate them to a large distance, we obtain just two hyperons. However, if the 
currents hi — h^, are represented in the form "singletC^singlet" with respect to color, they will also 
have the flavor structure "singletCgisinglet" (each color singlet is also a flavor singlet), and therefore 
their structure seems less natural than the structure of hg. 

As a result, we conclude that the current h^ is more "physical" than the currents hi — h/i, since it has 
a nonrelativistic limit and is constructed as "octet^octet" with respect to flavor (in the sense indicated 
above). However, the currents hi — h^ are "poor" because they either do not possess a nonrelativistic 
limit (the currents hi, and /14) or have the flavor structure "singlet^singlet" (the currents hi,h2, 
and /13). The current /15 (we recall that the results for it are practically identical to the results for 
/le) has the flavor structure "nonetCgmonet" , and the "octet(8)octet" term in it dominates over the 
"singlet(Xisinglet" term. 

We must make the reservation that formally the method of QCD sum rules dictates a priori 
the only condition for the choice of the current: it must possess the required quantum numbers. 
However, a posteriori (after the fit) positive results were obtained only for the current satisfying the 
two requirements formulated above. In our opinion, this indicates the existence of a criterion which 
makes it possible to select the optimal (" physical" ) current from the set of currents with the quantum 
numbers of the given channel. We note that this conclusion is confirmed also by the results of Ref. |15j . 
in which the deuteron mass was calculated by the same method. In that paper, the current was also 
constructed as "baryon(8)baryon" , and positive results were obtained only when in the current the 
term having a nonrelativistic limit dominated over the term not having this limit. 

Unfortunately, the relatively low accuracy of our method in the determination of the mass (~ 15%) 
does not make it possible to draw a conclusion on the basis of our values (14,15) about whether the 
mass of the dihyperon H lies below or above the AA threshold (equal to 2.23 GeV). (We recall that 
the high accuracy of the MIT quark-bag model permitted Jaffe to conclude that H lies below the AA 
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threshold and is therefore stable.) On the basis of our analysis, we can say that the method of QCD 
sum rules agrees with the possible existence of a dihyperon H either below or above the AA threshold 
(since in both cases it is possible to represent the spectral density p in the form of the sum of a pole 
contribution from H and a contribution from an effective continuum). 

We shall consider in detail the case in which H lies below the AA threshold. In this cases, if the 
binding energy e = 2m\ — mn is sufficiently small (what " small" means will be particularized below) , 
there can occur the interesting effect of resonance enhancement of the contribution of the AA state 
of the continuum T^. To describe this effect, let us consider the contribution of the state |AA) (two 
free hyperons) to the spectral density: 

PAA(g') = {27Tf{0\h{0)\AA){AA\h+{0)\0)5{q-pi-p2) (19) 

where Pi = P2 = '^aj ^^'^ ™ integration over the phase space of the AA state is understood. (We 
note the inclusion of the contributions of the and ^.N states in the approximation of exact SU^{3) 
symmetry with allowance for only the flavor-singlet states does not alter the results.) For the matrix 
element (0|/i(0)|AA) = fh{s), s = q^, we can write the dispersion relation 

fhis) = ^^^^ + f ^ subtractions (20) 

'''^ ' m^-S Jirnl s' -s ^ ^ 

where Xh is defined in (2) and qha is the constant for the interaction of the dihyperon H with AA. 
On the other hand, for energies which are not very large, when only the elastic channel of scattering 
of A by A is possible (we do not take into account the weak interactions), unitarity gives 

h{s) = {l + 2ikTj,^{k))ms) (21) 

where Taa = ~l/('* + i^) is the amplitude for elastic resonance AA scattering, k = ^/fnXe (e = 
2mA — fnil), s = 4:(m? + k'^), and k = |k| is the momentum of the A hyperon in the c.m.s. Comparing 
the pole behaviors of Eqs. (I^Uj) and (PT|) . we obtain the expression for the pole of fh{s): 

Substituting (j^^ into (|T!?|) . we find the required expression for the contribution of the state |AA) to 
the spectral density: 

1 X^g^j^ 1 1 /s-4m|\ 

The interaction constant gnA can be found \16^ by comparing the pole contribution of the virtual 
dihyperon H to the quantum field-theoretical expression for the amplitude of elastic AA scattering 
with the pole behavior of the quantum-mechanical expression for the amplitude of low-energy elastic 
resonance AA scattering: 

g%j, = 327r4mAK . (24) 

It can be seen from Eq. (|23j) that paa('S) has a resonance peak near the AA threshold with a maximum 
at the point Smax = ^tti^ + 4k^ and a width determined by the value of k. Thus, for small k (when 
the width of the resonance peak does not exceed the error in the method of QCD sum rules, i.e. 15% 
of m}{) this peak and the peak from the dihyperon H itself fuse in the spectral density p into a single 
effective peak, the position of whose maximum we determine by fitting the sum rule. However, since 
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the width of this peak for small k does not exceed the error of the method in the determination of niH, 
by calculating the position of its maximum we are calculating, within the accuracy of the method, the 
value of rriH itself. However, the constant A in the spectral density (9) will then no longer be equal to 
the projection Xh of the current h onto the dihyperon state [see (2)]. This will be some effective A, 
whose relation of Xh we shall now estimate. 

For this, we estimate the contribution of the resonance peak from the AA state to the sum rule in 
the leading order in k. Using Eqs. (|23|) and (|24j) . we find 

r e-^/ypAA{s)ds = e-yyxl^-^^^ . (25) 

Here sq = (2mA + -E'o)^ must be taken such that the range of integration spans the resonance peak, 
i.e. Eq is approximately equal to the width of the peak. Thus, the presence of the resonance peak 
from the AA state leads to a difference between the effective constant A in Eq. (9) and the constant 
Xh for coupling of the current h to the dihyperon H] 

X^ = Xjj{l + Tr-^^mAEo/K). (26) 

When K is fairly large (i.e. when the width of the AA resonance peak is larger than the uncertainty 
in the method of QCD sum rules in the determination of the masses), the peak becomes quite weak 
and can be included in the effective continuum. 

Finally, we formulate our main conclusion. The analysis of our sum rules favors the existence of a 
dihyperon H, in other words, it agrees with Jaffe's prediction in the framework of the MIT quark-bag 
model that there exists a dihyperon H with the quantum numbers of two A hyperons. However, the 
accuracy of the method of QCD sum rules, which in our case is ~ 20%, does not make it possible to 
determine whether the mass of the dihyperon H lies above or below the AA threshold (i.e. whether 
H is stable). 

The authors are grateful to A. N. Tavkhelidze for interest in the work and to K. G. Chetyrkin for 
helpful discussions. 
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